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Microlensing Evidence for Super-Eddington Disc Accretion in 
Quasars 



(N 

o 

(N 
D 



O 

U 

Of 

6 
u 

> 
oo 
t> 

00 

o 

(N 



P. Abolmasov* and N. I. Shakura, 



Sternberg Astronomical Institute, Moscow State University, Universitetsky pr, 13, Moscow, 119992, Russia 



Accepted — . Received — ; in original form — 



1 INTRODUCTION 



Since the work of iLvnden- Bell {1969), disc accretion onto super- 
massive black holes is a commonly accepted interpretation for the 
activity of quasars, both radio-loud and radio-quiet (sometimes dis- 
tinguished as quasi-stellar objects, QSO). Among all the active 
galactic nuclei, quasars are distinguished by higher luminosities 
(exceeding that of host galaxies) that is most likely connected to 
higher accretion rates. 

Spectral energy distributions in optical and UV are rea- 
sonably consistent felvis et all 1 1994) with the predictions of 
the standard thi n accretio n disc m odel introduced in th e sem - 
inal works of IShakural Jl972h; [shakura & Sunyaevl < fl973h ; 
iLvnden-Bell & Pringlej dl974l) : iNovikov & Thornel Jl973h . How- 
ever, the predicted angular sizes of quasar accretion discs are too 
small (microarcseconds and less) to be resolved directly. For to- 
day, quasars remain essentially point-like ("quasi-stellar") objects 
resolved only indirectly, in particular by mic rolensing effects. 

As it was shown by lAgol & Krollkl Jl999h . microlensing 
by the stellar population of lensing galaxies is sensitive to the 
size of the emi t ting region. Here, we adopt the statement of 
iMortonson et alj J2005h that the basic quantity that microlensing 
amplification maps and curves are sensitive to is half-light radius. 
Half-light radius R 1 /2 is defined as the radius inside which half of 
the observed flux is emitted at a given wavelength. 

Numerous studies aimed on probing the spatial properties 
of quasar accreti on discs with h elp of microlensing. While most 
early works (see IWambsganssI l|2J)0g) for review) reported rea- 
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ABSTRACT 

Microlensing by the stellar population of lensing galaxies provides an important opportunity 
to spatially resolve the accretion disc structure in strongly lensed quasars. Disc sizes esti- 
mated this way are on average larger than the predictions of the standard Shakura-Sunyaev 
accretion disk model. Analysing the observational data on microlensing variability allows to 
suggest that some fraction of lensed quasars (primarily, smaller-mass objects) are accreting in 
super-Eddington regime. Super-Eddington accretion leads to formation of an optically-thick 
envelope scattering the radiation formed in the disc. This makes the apparent disc size larger 
and practically independent of wavelength. In the framework of our model, it is possible to 
make self-consistent estimates of mass accretion rates and black hole masses for the cases 
when both amplification-corrected fluxes and radii are available. 
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sonable agreement between the observational data and the stan- 
dard accretion disc theory, several important results are at odds 
with the theoret ical predictions. St udying microlensing amplifi- 
cation statistics, IPoolev et alj d2007h find best-fit disc sizes more 
than one order of magnitude larger than the theoretical predic- 
tions based on photometrical data. Partially, this may be at- 
tri buted to the mass estimates u sed in this study (see di scussion 
in I Abolmasov & Shakuri] j2012h ). In lMorgan et alj ( |2010|) . incon- 
sistency is somewhat smaller (about a factor of 3) but still sig- 
nificant. Accretion discs seem too large for their apparent lu- 
minosities or too faint for their sizes in the UV/optical range 
(~ 2000-4000A) Many pap ers (such as j jimenez- Vicente et alj 
d2012l) . |Poolev et al.l ( [200% and lMorgan et alj fcOld) ) interpret this 
inconsistency as an evidence for insufficiency of the standard ac- 
cretion disc model, but no universal solution was proposed so far 
to account for the size discrepancy. There are i ndications for possi- 
bly higher black hole masses in some objects ( iMorgan et aHl2010l : 
lAbolmasov & Shakurall2012l) , but no changes in masses, accretion 
rates and efficiencies can explain the observed sizes and fluxes si- 
multaneously. 



One of the important issues in quasar microlensin g studies 
is whether the disc radial scale Rs dMorgan et alj 120101) depen- 
dence on wavelength is consistent with the power law Rs oc A 4 / 3 
predicted by the standard accretion disc theory. Several important 
disc models predict power law dependences R(X) oc A^ with dif- 
ferent exponents. Below we will refer to £ as "structure parame- 
ter". While for QSO J2237+0305, classical C = 4 /3 works fairly 
good jEigenbrod et all 120081 : lAnguita et all l2008t) , other objects 
such as SPSS 092 4+0219, for instance, clearly require smaller f, 
lFlovdetalJ ( l2009l) propose angular momentum inflow at the inner 
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edge of the disc in SDSS 0924+0219 dAgol & Krolikl koOO) as an 
explanation for the apparently very steep temperature law in the 
disc. This model implies Rs oc X 8 ^ 7 marginally consistent with 
the observational data. 

In the recent work bv lBlackburne et al.1 JioTTh . several other 
objects were shown to have much shallower i?(A) dependences, 
some consistent with Rs = const for a broad range of comoving 
wavelengths, O.lfi < A < lfi. The only object having conventional 
thin-disc scaling is MG J0414+0534 that has the highes t mass 
among the sample of 12 objects considered by iBlackburne et aU 
d201 lh . All the smaller-mass (M < 10 9 M o ) black holes are char- 
acterised by Q ~ 0-0.5. 

Microlensing effects in the X-ray range are more profound 
than in the optical dPoolev et al]|2007h . Independently of the disc 
structure in the optical and UV ranges, X-ray properties are more 
or less similar for all the obje cts where microlensing effects were 
studied in the X-r ay range dChen et al.l 1201 ll ; iDai et alj |2010| : 
Morgan et al. 2012). Evidently, X-ray e mission comes fro m some- 
where inside the inner ~ 10 x GM/c 2 Jchen et alj20l3) . and the 
exact mechanisms driving the formation of the X-ray continuum 
and lines are yet to be revealed. 

Small structure parameters originate not only for very steep 
temperature slopes in multi-blackbody models. For instance, C = 
is naturally reproduced if the brightness distribution does not de- 
pend on wavelength. This may be achieved if the accretion disc is 
surrounded by an envelope optically thick to Thomson scattering. 
Without affecting its spectral properties, scattering changes the spa- 
tial brightness distribution of the disc radiation. In general, accre- 
tion disc will increase its apparent radius and lose its intrinsic radius 
dependence on wavelength. A possible origin for such a scattering 
envelope is super-Eddi ngton accretion that leads to formation of 
a Thomson-thick wind dShakura & Sunvaevll 19731) . Since there is 
observational evidence that some qu asars accrete in s upercritical 
regime, especially at larger redshifts dCollin et alj|2002l) , we con- 
sider this scenario quite plausible. 

In the following section[2] we describe a simple scattering en- 
velope model that we use to account for the spatial properties of 
microlensed quasars. It will be shown that such an envelope may re- 
sult from super-Eddington accretion by a supermassive black hole. 
In section[3]and[4] we describe the observational data and interpret 
them in the framework of the scattering envelope model. In section 
[5] we consider the possible connection between the putative class 
of supercritical quasars, broad absorption line (BAL) quasars, and 
make conclusions in section[6] 



2 SPHERICAL ENVELOPE MODEL 

Though the issue of super-Eddington accretion is complicated, 
and different effects like photon trapping should be taken into 
acc ount, the simple pic t ure of supercritical accretion introduced 
by IShakura & Sunvaevl i 19731) is sufficient for our needs. This 
picture is supported by comprehensive num erical simulations 
dOhsuga et al.ll2005l : IOhsuga & Mineshigell201 lh . 



2.1 Transition to supercritical regime 

Radiation pressure is the principal feedback source for disc ac- 
cretion. While for a spherically-symmetric source, both radiation 
pressure force and gravity scale oc R~ 2 with distance that leads 
to the universal Eddington luminosity limit, accretion disc geome- 
try makes the situation more complicated because the source is no 



longer isotropic, and the forces are no longer collinear. For a thin 
disc, vertical gravity component grows with the vertical coordinate 
z while the flux generated in the disc does not significantly depend 
on z high enough in the atmosphere. Thin accretion disc supported 
by radiation pressure will have thickness determined by the equi- 
librium of the vertical components of the two forces. 



x _ x 3 GMM 
~c~ ~ R 3 




GM 



H 



Here, R is the radial coordinate, Ri n is the inner disc edge that 
for a black hole disc it is instructive to identify with the innermost 
stable orbit radius, x ~ 0.4cm 2 g _1 is electron scattering opacity, 
c is the speed of light, G is gravitational constant, M and M are 
the black hole mass and accretion rate. H is accretion disc half- 
thickness that may be subsequently expressed as: 



3 xM 
8^~c~ 




3 GM 



M* = 



Here, we normalised the mass accretion rate as M — M*rh, 
L E dd 4ttGM 



(1) 



where Lsdd = 'inGMc/ x is the Eddington luminosity. 
When the thickness of the radiation-supported disc becomes com- 
parable to its radius, thin-disc approximation breaks down and the 
balance of forces is inevitably shifted: gravity scales as (z 2 + 
R 2 ) -1 , while the disc radiation decreases more slowly because 
at these distances accretion disc is still a strongly extended ra- 
diation source. Hence we assume the disc super-Eddington if its 
equilibrium half-thickness is H ^ R. This condition defines the 
spherisation radius in a non-relativistic regime but with a correc- 
tion term that may play non-negligible role near the critical accre- 
tion rate. Let us introduce dimensionless radius r = R/Ri n , where 
R GM/c 2 is the inner radius of the disc. Dimensionless 

inner radius Xi„ varies between 1 (extreme Kerr case, corotating 
disc) and 9 (extreme Kerr, counter-rotation). Existence of a super- 
critical region in the disc requires existence of a root r > 1 of the 
following equation: 

3 ™ (2) 



1 - 1 



2 Xi n 



This equation may be reduced to a cubic equation for <Jr hav- 
ing either two or no positive real roots, depending on the right-hand 
side. The minimum of the left hand side is always at r cr = 9/4. 
This implies the critical accretion rate of rh cr = i.5xi n . The lu- 
minosity of a non-relativistic disc with this value of mass accretion 
rate is L/LEdd = m cr r\ ~ 9/4, where r\ = L/Mc 2 is accre- 
tion efficiency. This value may be thought of as the gain that disc 
geometry provides for the Eddington limit. Without the correction 
term, the critical accretion rate should be about an order of magni- 
tude lower, m cr ,o = 2xi n /3. This justifies the attention we pay to 
the correction term. Still, the real transition to supercritical accre- 
tion may be much more complicated and influenced by relativistic 
effects and the exact physics of wind acceleration. 

Spherisation radius for m > m cr may be found in a way that 
takes into account the correction term: 



3 GM . 2 / • / \ 



(3) 



Here, ip is the correction multiplier that accounts for the influ- 
ence of the inner radius in the disc. It may be found as the largest 
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root of the cubic equation for y/r that follows from l|2}. When ac- 
cretion is supercritical (m > rh cr — 4.52;™), the cubic equation 
has three root s and its solutio n may be expressed using trigonomet- 
ric functions dNickallsll 19931) . Solving the cubic equation yield the 
following expression for correction multiplier: 



^ )= ^ cos G arccos ("^ 

For m 2> 1, ip(rh/xi n ) rapidly approaches 1, hence: 



(4) 



R 



sph 



3 GM 
2 mX — ' 



that coincides with the cla ssical definition of spherisation ra- 
dius jShakura & Sunvaev|[l973h . 

Formation of disc winds in supercriti cal accretion regime 
is supported by numerical simulations dOhsuga et al] 2005; 
Okuda et alj|2005h . If the radiation flux from the disc exceeds the 
local Eddington value, some part of its energy is converted to ki- 
netic energy of the outflow. Terminal wind velocity v w may be es- 
timated through Bernoulli integral: 



GM 
~R~ 



- 1 



For example, setting the luminosity to the "disc Eddington 
limit" L = jLEdd implies v w = ^2.5GM/ R m where R w is 
some effective radius where the outflow is formed. Since most of 
the matter is ejected from R ~ R sp h, we assume the terminal 
velocity of the wind is proportional to the escape velocity at the 
spherisation radius: 



v w = P^2GM/R aph , 



(5) 



were /3 ~ 1 is some additional dimensionless multiplier intro- 
duced to account for the uncertain details of wind acceleration and 
formation of the outflow. 



2.2 Spherical envelope radius 

Let us assume that the envelope is composed of a fully ionised 
spherically-symmetric wind expanding at a constant velocity. Con- 
tinuity equation allows to connect electron density with the mass 
loss rate M w = f w rhM* and outflow velocity v w calculated ac- 
cording to equation (|5j. 



P = 



<inR 2 v n 



Envelope size is defined by the radius where the radial optical 
depth toward the observer is unity. 



r(R) 
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xp(R)dR ■- 



fwxM 
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envelope has a radius practically independent of wavelength while 
the spectral properties of the scattered radiation remain more or less 
unchanged. The envelope is actually a pseudo-photosphere: it ex- 
pands supersonic ally at a large, possibly mildly relativistic velocity 



2.3 Apparent intensity distribution 

Half-light radius is defined by the general relation that may be used 
for any radially-symmetric intensity distribution: 



J* 1 / 2 I(R)RdR _ i 
J+°° I(R)RdR = 2 



(7) 



For a standard thin non-relativistic accretion disc, monochro- 
matic intensity scales with radius as: 



1 



exp (R/R s ) 



3/4 



i(R) 



where Rs is the disc radial scale defined by condition hv £m = 
kT(Rs) without the correction factor. If R in <C R1/2, to a high 
accuracy: 



R (disc) _ 2 .URs 



~ 2.44 
= 2.44 



45X* m GMM 
16w B hc' 2 
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(8) 



4n b h>eGM 
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Here, A em is the comoving (quasar reference frame) wave- 
length of the observed radiation {v em = c/A em is corresponding 
frequency; observed wavelengths and frequencies are denoted as 
Xobs and Vobs), h is Planck constant. 

For a spherical envelope, brightness is nearly uniform in the 
centre and declines as a power law at large radii. We have taken the 
extended scattering atmosphere model described in Appendix [A] 
and calculated the intensity at infinity for different shooting pa- 
rameter values coming to the overall conclusion that the half-light 
radius in this model is proportional to the photosphere radius and 



R 



{envelope.) 

1/2 



l.O&Ri. 



2.4 Disc radiation 

Standard accretion disc temperature law may be written (neglecting 
the correction term important for the inner parts of the disc) as: 



T { R)=(I^R-^ 
1 2 axe 



1/4 



where a is Stephan-Boltzmann constant. Monochromatic flux 
is found as an integral over the picture plane: 



Ri = R(t = 1) = 



fw*cM 
4ttv w 



R\c 2 j3 f w . 3/2 , , ■ 1 
~GM ~ V 8~/3~ m 



(6) 



For a moderately super-Eddington accretor with f w ~ 1, the 
size of the envelope becomes comparable and may even exceed the 
size of the accretion disc in the ultraviolet range. Such a scattering 



F u 



I v dQ, 



2tt cos i 



D 2 x (1 + z) 



I v {R)RdR 



We use angular size distance D — D(z), and i is disc incli- 
nation. Assuming the radiation generated in the disc has locally a 
blackbody spectrum leads to the following estimate for monochro- 
matic flux valid far away from the high- and low energy cut-offs 
connected to the inner and the outer disc edges: 
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F v = 8tt 



'2\l/3 



r(8/3)C(8/3) 



(GM) 4 



c 8/3 h 5/3 K 2/3 CT 2/3 ' 

1 

D 2 X(l + z) 8 / 3 



(9) 



Here, T and C, are Gamma-function and Riemann zeta- 
function (not to be confused with the structure parameter £ that 
is used without any argument). The above formula may be used to 
estimate the mass of the supermassive black hole (SMBH) using 
one observable quantity (flux) and one unknown parameter m. 



M 



x (87rr(8/3)C(8/3))- 3/4 x 
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Gk^lii (io) 
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Spherical envelope scrambles the radiation generated in the 
disc and makes it roughly isotropic. Effective inclination cosine 
in this case is cosi e // = 1/2 since the initially anisotropic flux 
F v oc cos i is re-distributed isotropically with the total luminosity 
conserved. 



2.5 Mass and accretion rate estimates 

Microlensing studies are unique for distant quasars since they al- 
low to estimate the size of the emitting region in continuum inde- 
pendently of the observed flux. In the case of a standard accretion 
disc both observables, F v and R1/2, may be used to estimate only 
one combination o f black hole mass and a c cretio n rate, M 2 rn (see 
also discussion in lAbolmasov & Shakura! fe012h ). Existence of a 
scattering envelope allows to break this degeneracy and make self- 
consistent estimates of both principal parameters (M and m). 

Solving the system of two equations d 1 0| > and ® for M and 
m allows to estimate both black hole mass and dimensionless ac- 
cretion rate. We also set cosi = cosieff = 1/2: 



rh = (87rr(8/3)C(8/3)) 3/4 
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Here, we expressed the observed flux through the magni- 
tude I in the HST F814W band (\ obs ~ 0.79/Lt) adopting the 
zero-point flux equal to F v = 2.475 x lCT 20 ergcm~ 2 s" 1 Hz -1 



(Holtzman et al. 995 ) since we use the amplification-corrected 
magnitudes from Morgan et al.l {20105 obtained in this photomet- 
rical band at the HST. 

Since the left-hand side of dl lb dependence on the mass ac- 
cretion rate is much stronger, simple forward iteration works good. 
Once m is found, the black hole mass may be estimated following 
flOll as: 



R 



(disc) 

1/2 



R 



(envelope) 

1/2 



2.44i? s = 1.06i?i 

After substituting equations {8} and ([6j, one gets the following 
mass limit: 



MUm = (M) 3 (|) 3/2 # (£) 3 X 

~3.8x 10 10 (^-) 4 (^)- 7/2 ^- 3 (m/ a;i „)M 



(12) 



For higher masses (at a given wavelength, for fixed m), the 
size of the envelope is smaller than the half-light radius of the 
disc, and the appearance of the quasar will be close to the thin 
disc case. The limit is shown in figures [T] and [2] with solid lines. 
This limit ev i dently depends on wavelength. For the sample of 
Mor gan et al] ^2010), comoving frame wavelength changes in the 
range 0.2-0.5/i that corresponds to about a factor of 2 upward 
shift of the limit in figureQ] 



3 MULTI-WAVELENGTH DATA 

3.1 Disc radii and amplification-corrected fluxes 

We make use of the amplification-corrected fluxes and 
microl ensing-based radii collected and publish ed bv lMorgan et al.l 
fcOlOh. The sample of iMorgan et alj J2010h overlaps with this 
of iBlackburne etaL In the original work, disc radii are 

given in terms of accretion disc radial scale Rs- Since we propose 
that at least for some objects, emitting region has different nature 
and intensity distribution, we recalculate these disc radii into 
model-independent half-light radii. Fitted with a standard-disc 
model with radial scale Rs as defined above, emitting region may 
be characterised by the half-light radius of R1/2 — 2.44Rs- The 
scattering envelope model we use has one characteristic radius Ri 
where optical depth equals unity. For this model, R1/2 — I.O6-R1 
(see Appendix lAt . Ha lf-light radii are give n in table 14. 1| For all 
the ob jects studied by Morgan et al. (2010) and Blackburn e et al.l 
J201 ll) . the two estimates are consistent within the uncertainties. 
Since the uncertainties are very high, it is difficult to set any 
constraints upon the possible var iations of the half-lig ht radii. 

For their sample of objects. IMorgan et alj ( 120101) also provide 
magnitudes corrected for strong lensing amplification (see also ta- 
ble l4.lt based upon HST observ ations in the F814W fil ter. Flux cal- 
ibration is based on the paper of lHoltzman et al.l (Tl995). see section 



M ~ 4.6 x 10 



7 x 10 -o.3(/-i9) ( £(zl\ 3/2 + ^2^-1/2 

\iGpcy 



Mr: 



2.6 Disc and envelope sizes 

Depending on the wavelength range, a supercritical disc may be ei- 
ther observed directly (if its size is larger than the photosphere of 
the wind) or covered by the photosphere of the supercritical wind. 
Equality of the half-light radii set by the two radial scales leads to 
the following condition for observational importance of the enve- 
lope: 



3.2 Masses and emissivity slopes 

Several methods are used to estimate masses of supermassive black 
holes. Most of them are model-dependent and suffer from biases of 
different nature. For bright distant quasars, masses are usually esti- 
mated either through photometrical data (bolometric luminosity is 
restored from multi-wavelength observations) or by measuring the 
widths of broad emission l ines and the size of the emi tting region 
by reverberation mapping (Blan dford & McKed 1982). While the 
first method relies heavily on ad hoc assumptions about the mass 
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Table 1. Properties of the microlensed quasars from the sample of Morgan et al, 12010). Black hole mass Mg and accretion rate rn are calculated for the case 
of supercritical accretion for a = 0.9. 



Object 



M^jIO^Mq I corr ,mag 



R 1/2 , 10 15 cm 
(Morg an et al. 2010) 



R 1/2 , 10 15 cm M S ,10 9 M Q 
(Blackburne et al. 2011) 



HE J0435-1223 0.50 20.76±0.25 

SDSS 0924+0219 0.11 21.24±0.25 

FBJ0951+2635 0.89 17.16±0.11 

SDSS J1004+4 112 0.39 20.97±0.44 

HEJ1104-1805 2.37 18.17±0.31 

PG 1115+080 1.23 19.52±0.27 

RXJ 1131-1231 0.06 20.73±0.11 

SDSS Jl 138+0314 0.04 21.97±0.19 

SBSJ1520+530 0.88 18.92±0.13 

QSO J2237+0305 0.9 17.90±0.44 

QJ0158-4325 0.16 19.09±0.12 

O J0158-4325 iMorgan et alj2012h 0.16 19.09±0.12 
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accretion rate and accretion efficiency, the second has a fundamen- 
tal uncertainty connected to the geometry of the emitting region. 
Virial mass is estimated as: 



M = f- 



'Rblr 
~G ' 



where a 2 is the velocity dispersion corresponding to the ob- 
served line width, Rblr is the size of the emitting region (deter- 
mined with help of reverberation mapping), and coefficient / is cal- 
ibrated using better -studied nearby active galaxies where / ~ 5.5 
dOnken et al]|2004h . Sometimes only a limited number of spec- 
tra is available and reverberati on analysis in impossible. In this 
case, empirical virial relations iVesterg aard & Peterson 20(3) ar e 
used. These two types of mass estimates will be hereafter referred 
to as virial. Since w e use the microlensing-based disc radii from 
iMorgan et sul fa) 10) we also make use of the virial m asses given 
in this work (see references in this paper, especially IPeng et all 

Jiooi). 



4 RESULTS 

4.1 Masses and accretion rates 

Mass and dimensionless mass accretion rates estimated from the 
observables by the method introduced in section [231 are given in 
table |4. l| (for a = 0.9 and Xi n = 2.32) and shown in figures[T]and|2] 
for two values of a and Xi n . All the masses and mass accretion rates 
were determined in the assumption of existence of an optically- 
thick scattering envelope. They apply only to the objects where the 
disc is surrounded by an envelope larger than the disc itself (i. e. 
above both the dotted and the solid lines in figure[]}. 

Errors given in table and figures[T]and[2]were calculated using 
direct n on-linear error prop agation. We used the lcr uncertainties 
given bv lMorgan etal] hold) for the radii and fluxes. Since we do 
not know the exact probability distributions of these quantities, it 
seems to be the most reasonable approach. We substituted Ri ± 
ARi and I ± Al into equations dl lb and J 1 0b and interpreted the 
highest and the lowest values of m and M as the ends of some 
representative confidence interval. 

Properties of the larger part of the objects shown in figure[T]are 
consistent with accretion in a moderately super-Eddington regime. 



If the accretion efficiency is high (a > 0.9), all the objects may 
be interpreted as super-Eddington. Most of the objects are how- 
ever difficult to identify as super- or sub-Eddington sources due to 
large uncertainties in rh. More probable super-Eddington objects 
such as RXI 1131-1231 and PG 1115+080 tend also to have lower 
( (see section l4~2l . Einstein's cross, one of the most probable sub- 
critical discs from the sample, is evidently among the lowest-m ob- 
jects. According to our model, Q J0158-4325 should best conform 
to the thin dis c model. Indeed, disc size for this object reported in 
IMorgan et ailbOlCl) is in good agree ment with the theoret ical pre- 
dictions. However, the recent work of IMorgan et alj reports 
evidence for a larger disc size, several times larger than (but still 
consistent at about 1.5<r confidence level) the standard model pre- 
dictions based on the measured amplification-corrected flux. Since 
dimensionless mass accretion rate is proportional to the envelope 
radius, the larger radius makes the properties of Q 10158-4325 con- 
sistent with our supercritical disc model. 

Masses determined in the spherical envelope model are gener- 
ally smaller than virial masses (see figure^. Indeed, applicability 
of the classical virial relations to an expanding supercritical wind is 
questionable. As we will show in the section[5] the broad emission 
lines observed in quasar spectra are unlikely to be formed in the 
outflowing matter, but one still may expect violations of the virial 
relations derived for sub-critical active galactic nuclei. 



4.2 Structure parameter correlation with mass 

Several quasar microlensing studies used multi-wavelength data to 
trace the disc size dependence on wavelength. Fitting this depen- 
dence with a power law allows to check the validity of several ac- 
cretion disc models such as thin disc (£ = 4/3), slim or irradiated 
disc (£ = 2) and a thin disc with a strong torque at the inner radius 
(£ = 8/7). Available data are collected in tableland in figure|4] 

Interestingly enough, all the objects where standard disc slope 
works have masses M > 1O 9 M . All the smaller -mass black holes 
show intensity distributions with much lower £. Besides the large 
fitting errors, objects in figure [4] may be separated in two groups: 
high-mass black holes surrounded by accretion discs similar to 
standard and lower-mass objects where £ ~ 0-0.5. An evident 
qualitative solution is to propose that at least some quasars accrete 
in super-Eddington regime. If all the bright lensed quasars accrete 
at M ~ 3OM0 yr _1 , Eddington luminosity will be reached for 
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M ' O g k 



0.10 
M s , 10' M soklr 



Figure 1. Quasar parameters (mass and dimensionless mass accretion rate) recovered in the framework of supercritical envelope. The left and the right 
panels correspond to a = (Xi n = 6) and a = 0.9 (Xi n ~ 2.32), respectively. Horizontal dotted line shows the actual Eddington limit: accretion is 
sub-critical below the line and super-critical above it. Solid inclined line marks the limit whe re the envelope becomes larger than the disc monochromatic size 
(for A e m = 0.25/j) . Hereafter, thick dots show objects that are also present in the sample of Blackburne et al. (2011), and the recent result for Q J0158-4325 
iMorgan et al. 2012) is shown by a star. Errors correspond to ltr uncertainties in flux and radius. 





M s , 10' M sotar M s , 10' M so i ar 

Figure 2. Same as previous figure, but instead of dimensionless mass accretion rate fa, dimensional mass accretion rate M is given. 



M ~ 1O 9 M0 for accretion efficiency rj ~ 0.1. Lower-mass ob- 
jects are expected to enter super-E ddington accretion phase e asier 
and lose excess accreting matter dShakura & Sunvaevlll973l) . We 
propose that supercritical wind does not affect the spectral energy 
distributions of QSO much but changes its spatial properties, acting 
as a lampshade that changes only the visible size and shape of the 
lamp. 



v ~ m -1 / 2 c. At the same time, Doppler widths of broad emission 
lines are significantly lower, v/c ~ 0.01. Besides, the expected 
emission line luminosities are several orders of magnitude less than 
that of observed broad emission lines in quasars. If we propose a 
constant filling factor /, a recombination line in a supercritical out- 
flow should have a luminosity estimated as the following volume 
integral: 



5 DISCUSSION 

In previous sections we have shown that some lensed quasars 
are surrounded by scattering envelopes of moderate optical depths 



Rx/R sph 



,1/2 



1 - 10. Expected outflow velocities are 



Lune = —ahv x / n e mdV 
J Jv 

Here, a is the recombination coefficient, n e and n, are the 
electron concentration and concentration of the particular ion emit- 
ting the line. It is convenient to express the concentrations as 
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m — n e Xi = pxi/m, Below, we fix the values of Xi and the 
effective particle mass m. For completely ionised hydrogen-rich 
gas, m is about the proton mass. 

L Une ^ w -7 fw^ifn 2 OL fOOOA R sph 

L E dd ~ f/3 2 10 _13 Cm 3 S _1 Xline Rin 

Integration is performed from some inner radius Ri n to in- 
finity. Since the inner parts of the flow are considerably ionised 
the above luminosity is an upper estimate. Predicted equivalent 
widths are of the order ~ 1CP 4 A, about five-six orders lower than 
the observed equivalent widths of broad emission lines. Condi- 
tions are much more favourable for formation of absorption lines, 
since the wind is thick to electron scattering and Nh ~ t jm v ~ 
lO^m^cm- 2 . 

Outflows are expected to manifest themselves in blueshifted 




Figure 3. Virial masses plotted versus mass estimates in supercritical wind 
assumption . Kerr parameter was set to a = 0.9. 

Table 2. Struct ure parameter as a func tion of SMBH mass. All the data 
were taken from Blackburne et al. 12011), if not stated otherwise. 



Object 


1O 9 M 


C 


reference 


MG J0414+0534 


1.82 


1.5 ±0.8 




MG J0414+0534 


1.82 


1.5 ±0.5 


Bate et al. (20081 


HE J0435-1223 


0.50 


0.7 ±0.6 


RXJ 091 1+0551 


0.80 


0.17 ±0.4 




SDSS 0924+0219 


0.11 


0.17 ±0.5 




SDSS 0924+0219 


0.11 


0.7 ±0.4 


Flovd et al. (2009) 


HEJ1104-1805 


2.37 


1.65 ±0.5 


Poindexter et al. (2008) 


PG 1115+080 


1.05 


0.4 ±0.5 




RXJ 1131-1231 


0.06 


0.4 ±0.5 




SDSS J1138+0314 


0.04 


0.4 ±0.5 




QSO J2237+0305 


0.90 


1.2 ±0.3 


Eieenbrod et al. (2008) 


QSO J2237+0305 


0.90 


1.1 ±0.3 


Eieenbrod et al. (2008) 






1 2+ 2 


(no velocity prior) 


QSO J2237+0305 


0.90 


Anguita et al. (2008) 



absorptions and P Cyg lines. Broad absorption line (BAL) quasars 
dTurnsheklll984T) show strongly blueshifted absorption components 
of UV and sometimes X-ray spectral lines. Two of the objects of 
our sample, PG 1 1 15+080 and SBS J1520+530, are known as BAL 
quasars. PG 1 1 15+080 also demo nstrates signatures of moderately- 
relativistic outflows. In particular, Chartas etHll d2003l) find two rel- 
ativistic absorption components (at ~ 0.1 and ~ 0.3c) of highly 
ionised iron species and an OVI absorption component at ~ 0.02c. 
Similar X-ray absorption systems were found for some other high- 
redsh ift luminous quas ars like APM 08279 +5255 JChartas et at] 
l2002h. HE J1414+1 17 dChartas et al]|2007l) and HS 1700+6416 
(Lanzuisi et alj|2012h . Absorption components with relativistic ve- 
locities were also found in the UV spectra of so me BAL, "mini- 
BAL" and narrow-absorption line quasars (see iNaravanan et al.l 
(2004) and references therein). Relativistic outflows often coexist 
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with slower absorption systems, and UV and X-ray absorption lines 
usually show different profiles and velocities. 

Discrepancy in wind velocities for different absorption com- 
ponents suggests that the wind is highly inhomogeneous, with dif- 
ferent components having different velocities. Besides, its structure 
may be far from spherical symmetry, and even in the spherically- 
symmetric case the shap e of the visible photosp here is distorted 
if the wind is relativistic dAbramowicz et alii 199 lh . For relativistic 
winds, there are two effects important for their spatial properties: 
(i) firstly, relativistic beaming makes the visible size of the emit- 
ting region ~ 7 3 times smaller, where 7 is Lorentz-factor, and ( ii) 
secondly, the optical depth along the wind flow direction is smaller. 
Both effects are expected to produce a wavelength-dependent limb- 
darkening effect that may be responsible for deviations of £ from 
for some objects. 

For several objects, the size of the X -ray emitting regio n 
was studied using microl e nsing effects ( s ee IPoplev et aO d2007t) : 
IChen et alj j201 ll |2012|) ; iMorgan et all d2012l) and references 
therein). Independently of the UV/optical structure parameter (, 
the X-ray sizes of all the studied quasars are estimated as several 
gravitational radii, that is sometimes one-two orders of magnitude 
smaller than the proposed envelope size. This is difficult to account 
for in the spherical envelope model since the size of the envelope 
(as well as the size of the accretion disc in the optical/UV range) is 
generally much larger. 

However, the outflows formed by super-Eddington accretion 
discs are expected to possess high net angular momentum that 
leads to for mation of an avoidance co ne also known as supercrit- 
ical funnel dShakura & Sunvaev|[l973l) . This picture is supported 
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Figure 4. Structure parameters for quasars of different masses. Virial 
masses are shown by crosses. For several objects(HE J0435-1223, 
SDSS 0924+0219, HE Jl 104-1805, PG 1115+080, RXJ 1131-1231, 
SDSS J1138+0314, QSO J2237+0305), we use our envelope model to esti- 
mate masses (shown by stars). 



by numerical simulations dOhsuga et al]|2005l ; lOkuda et alj|2005l) : 
for a large range of inclinations, the angular size of the source in 
X-rays is considerably smaller than the size of the outer photo- 
sphere of the wind. Optical depth of an envelope with a funnel 
is expected to be much lower (if the disc is viewed at low incli- 
nations, see IPoutanen et al.l d2007l) ) for the innermost parts of the 
disc where the X-ray component is supposed to be formed. The 
non-monotonic optical depth dependence on radius predicted for 
supercritical accretion disc winds can qualitatively explain the de- 
crease in the disc size at smaller wavelengths observed for some 
objects with C ~ suc h as PG 1115+080 and WF1 J2033-4723 
dBlackburne et al J201 il) . 

As long as we use virial masses, properties of th e X-ray radia- 
tion in dicate that it is formed near the last stable orbit dMorgan et al.l 
l2012h . However, smaller masses recovered in the framework of our 
supercritical envelope model shift the last stable orbit toward lower 
sizes. For instance, the X-ray emitting regi on of RXJ 1131-1231 
has the size of ~ 2 x 10 cm dPaietal.112 010). For the virial mass 
estimate of 6 x 1O 7 M0, this corresponds to ~ 7 x GM/c 2 , while 
the envelope-based mass is about an order of magnitude larger, and 
the estimated X-ray size becomes several tens of GM/c 2 . This dif- 
ference can hardly be used to distinguish between the individual 
mass estimates or individual models of X-ray emission separately. 
However, it should be borne in mind that virial masses are con- 
sistent with the models where X-ray emission is formed near the 
last stable orbit while the envelope-based mass estimates allow the 
X-ray emission to be extended for tens of gravit ational rad i i. This 
is consistent, for example, with the models like iLiu et alj d2012h 
where the X-rays are produced by the hot gas of the corona present 
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in the inner parts of the accretion flow. In any case, the size of the 
X-ray emitting region is expected to be considerably smaller than 
both the disc size at ~ 2500A and the scattering wind photosphere 
size. 

Existence of true absorption processes should also affect the 
apparent size of the spherical envelo pe. Assuming t otal th ermalisa- 
tion in a spherical relativistic wind, Fukue & lino] ( l2010h find that 
the visible photosphere surface follows the law T(R) oc R~ that 
implies £ ~ 1. More elaborate studies taking into account the tem- 
perature and ionisation structure of the wind are needed to explain 
the observed £ ~ 0-0.5 values of most of the putative supercritical 
quasars. 

True absorption processes are also important for wind acceler- 
ation. In the super-Eddington regime, wind is efficiently launched 
by resonance lines even in presence of strong X-ray radiation 
(Proga et al. 2000l; |Proga & Kallmanl2004T) . Resonance lines do not 
alter the measured photosphere size significantly since the wind is 
opaque to absorption only in a narrow wavelength range. How- 
ever, their contribution to wind acceleration through j3 and f w may 
be important. 

It is tempting to compare the population of super-Eddington 
quasars to the few known and well-studie d supercritical 
black hole X-ray binaries, primarily to SS433 iFabrikal 20041; 



Cherepashc huk. A. M. and Sunvaev. R. A. a nd Fabrika ~S. N. et al] 



20051) . For SS433, dimensionless mass accretion rate is of the 
order several thousands that implies much slower outflow velocity 
of ~ 1000km s -1 and thermalised emission from the wind 
pseudo-photosphere. On the other hand, mass accretion rates 
estimated in the p resent work, as well as the values given by 
ICollin et al.l d2002h . are considerably smaller. Maximal values 
are of the order m ~ 100-200. Note that these values are only 
moderately supercritical, Eddington luminosity is exceeded a fac- 
tor of ~ 10 - 20 (depending on the unknown accretion efficiency 
r\ ~ 0.06-0.4). It is more instructive to compare the population 
of super-Eddington quasars to the high-luminosity states o f X-ray 
binaries like GRS J1915+105 l lVierdavanti et al.l l2010h rather 
than to persistent stro ngly supercritical accre tors like SS433 or 
to sources like V4641 dRevnivtsev et al.l 120021) suffering strongly 
super-Eddington outbursts. 



6 CONCLUSIONS 

Scattering envelope formed by a super-Eddington accretion disc is 
a plausible model for the spatial properties of the emitting regions 
in some lensed quasars. Large spatial sizes (R ~ 10 16 - 10 17 cm) 
practically independent on wavelength are an expected outcome of 
a moderately super-Eddington (rh ~ 10 - 100) mass accretion rate. 
Black hole masses and mass accretion rates may be determined 
self-consistently if both disc size and flux estimates are present. 

The small sizes of X-ray emitting regions of microlensed 
quasars may be explained by existence of an avoidance cone, or 
supercritical funnel, in the disc wind. 

Some of our super-Eddington objects are broad absorption line 
quasars, and at least one (PG 1115+080) shows signatures of a 
mildly relativistic outflow. 
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APPENDIX A: BRIGHTNESS PROFILE MODEL 

Al Radiation transfer in Eddington approximation 

Optically-thick wind forms an extended pseudo-photosphere with 
a power-law density slop e, p oc R~ 2 . In general, radiation transfer 
is described by equation (Mihalas 197^): 



"55 M + 



i -ii 2 d 



{i) = -MR){s-i) 



R dfj, 

Here, I — I(R,fi) is monochromatic intensity, n = cos# 
is the cosine of the angle between the radius vector and radiation 
propagation direction, S = S(R) is source function. We consider 
only isotropic coherent scattering that allows to equate source func- 
tion to intensity averaged over solid angle. We use moment ap- 
proach, and use the first three radiation intensity moments defined 
as: 



■'-\ 



H = - 

2 



Id/j, 



Ifidfj, 



(Al) 



(A2) 



(A3) 



First radiation moment J has the physical meaning of mean in- 
tensity, hence S = J in our approximation. The system of moment 
equations is closed by Eddington's assumption K — fj where 
/ = 1/3 valid in diffusion approximation. Questionability of Ed- 
dington approxi mation for extended atmospheres is well known 
Chapm^ Jl966h but it is still sufficient for our purposes. Here, we 
consider pure electron scattering by a spherical atmosphere with 
electron density n oc p oc R- 2 . In this case, the system of moment 
equations takes the form: 



(A4) 



The system is simplified if / = 1/3 (inner parts, r 2> 1) 
and if / = 1 (opposite limit, r — > 0). Mean intensity (and hence 
source function) scales in these two approximations as oc CiR~ 3 
and oc C2R 2 X (1 + r), respectively. It is convenient to use the 
second formula and to set Ci — C2 = Ho = H(t = 1). Both 
asymptotics are then naturally reproduced. Ho may be connected 
to the physical flux at the photosphere as F(t = 1) = 4irHo 
and luminosity as L = (4tt) 2 R 2 Ho. However, this approximate 
formula does not conserve the total flux (intensity integrated over 
the solid angle deviates from the total radiation flux calculated as 
F(t = 1) = 4:-kHo) that may result in systematic errors, hence we 
adopt the following form for the source function: 



S{r) = Hot' 2 x (l + dx r~ 1/2 + r 



where r = R/Ri — 1/r, and d is a free parameter. Inte- 
grating source function for some shooting parameter P yields the 
observed intensity: 
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> 



Figure Al. Principal scheme illustrating integration along the line of sight performed in the Appendix. 



I = x J +0 ° S (\/P 2 +Z 2 ) e~ T(p ' l) p (Vp 2 + Z 2 ) dl, 
where t(P,1) is the optical depth along the current line of 



sight: 



T = xSl xP (VW+P~)dl = 

= (atan £ + tt/2) 

The coordinates and designations are shown in figure [ATI 
Finally, intensity distribution may be expressed as the follow- 
ing definite integral: 

I(p) = H x (w 2 (p) + u 3 (p) + dx u 6 /a(p)) (A5) 
where x — l/Ri andp = P/Ri, and: 



Jo 



-e/p ■ k 
e sin i 



Constant d is tuned in a way to fit the integral flux value, 
2n f + °° I(p)pdp = in Ho- Numerical integration allows to es- 
timate the value of d ~ —0.097. 

Half-light radius for this model may be estimated as -R1/2 = 
1.063.Ri to an accuracy of about 10 . Setting d = is also a 
reasonable approximation: it overestimates the flux by only about 
5%, while R1/2 — 1.05-Ri in this assumption. 



REFERENCES 

Abolmasov P., Shakura N. I., 2012, MNRAS, 423, 676 
Abramowicz M. A., Novikov I, D., Paczynski B., 1991, ApJ, 369, 
175 

Agol E., Krolik J., 1999, ApJ, 524, 49 
Agol E., Krolik J. H., 2000, ApJ, 528, 161 

Anguita T., Schmidt R. W., Turner E. L., Wambsganss J., Webster 
R. L., Loomis K. A., Long D., McMillan R., 2008, A&A, 480, 
327 

Bate N. E, Floyd D. J. E„ Webster R. L., Wyithe J. S. B., 2008, 

MNRAS, 391, 1955 
Blackburne J. A., Pooley D., Rappaport S., Schechter P. L., 2011, 

ApJ, 729, 34 

Blandford R. D., McKee C. E, 1982, ApJ, 255, 419 
Chapman R. D., 1966, ApJ, 143, 61 



Chartas G., Brandt W. N., Gallagher S. C, 2003, ApJ, 595, 85 
Chartas G., Brandt W. N., Gallagher S. C, Garmire G. P., 2002, 
ApJ, 579, 169 

Chartas G, Eracleous M., Dai X., Agol E., Gallagher S., 2007, 
ApJ, 661, 678 

Chen B., Dai X., Kochanek C. S., Chartas G., Blackburne J. A., 

Kozlowski S., 2011, ApJL, 740, L34 
Chen B., Dai X., Kochanek C. S., Chartas G., Blackburne J. A., 

Morgan C. W, 2012, ApJ, 755, 24 
Cherepashchuk, A. M. and Sunyaev, R. A. and Fabrika, S. N. et 

al. 2005, A&A, 437, 561 
Collin S., Boisson C, Mouchet M., Dumont A.-M., Coupe S., Por- 

quet D., Rokaki E., 2002, A&A, 388, 771 
Dai X., Kochanek C. S., Chartas G., Kozlowski S., Morgan C. W., 

Garmire G., Agol E., 2010, ApJ, 709, 278 
Eigenbrod A., Courbin E, Meylan G., Agol E., Anguita T, 

Schmidt R. W„ Wambsganss J„ 2008, A&A, 490, 933 
Elvis M., Wilkes B. J., McDowell J. C, Green R. E, Bechtold J., 

Willner S. P., Oey M. S., Polomski E., Cutri R„ 1994, ApJSS, 

95, 1 

Fabrika S., 2004, Astrophysics and Space 
Physics Reviews, v.12, p.1-152 (2004) 

(http://www.cambridgescientificpublishers.com/), 12, 1 

Floyd D. J. E., Bate N. E, Webster R. L., 2009, MNRAS, 398, 233 

Fukue J., lino E., 2010, PASJ, 62, 1399 

Holtzman J. A., Burrows C. J., Casertano S., Hester J. J., Trauger 
J. T, Watson A. M., Worthey G, 1995, PASP, 107, 1065 

Jimenez- Vicente J., Mediavilla E., Munoz J. A., Kochanek C. S., 
2012, ApJ, 751, 106 

Lanzuisi G., Giustini M., Cappi M., Dadina M., Malaguti G., Vi- 
gnali C„ Chartas G, 2012, A&A, 544, A2 

Liu J. Y„ Liu B. E, Qiao E. L., Mineshige S., 2012, ApJ, 754, 81 

Lynden-Bell D., 1969, Nature, 223, 690 

Lynden-Bell D., Pringle J. E„ 1974, MNRAS, 168, 603 

Mihalas D., 1978, Stellar atmospheres /2nd edition/. San Fran- 
cisco, W. H. Freeman and Co., 1978. 650 p. 

Morgan C. W., Hainline L. J., Chen B., Tewes M., Kochanek C. S., 
Dai X., Kozlowski S., Blackburne J. A., Mosquera A. M., Char- 
tas G., Courbin E, Meylan G, 2012, ApJ, 756, 52 

Morgan C. W., Kochanek C. S., Morgan N. D., Falco E. E., 2010, 
ApJ, 712, 1129 

Mortonson M. J., Schechter P. L., Wambsganss J., 2005, ApJ, 628, 
594 

Narayanan D., Hamann E, Barlow T, Burbidge E. M., Cohen 



© 2009 RAS, MNRAS OOO.rTHJJ 



Microlensing Evidence for Super-Eddington Accretion 1 1 



R. D., Junkkarinen V., Lyons R., 2004, ApJ, 601, 715 
Nickalls R. W. D., 1993, Mathematical Gazette, 77, 354 
Novikov I. D., Thome K. S., 1973, in Black Holes (Les Astres 

Occlus) Astrophysics of black holes.. Gordon and Breach, Paris, 

pp 343-450 
Ohsuga K., Mineshige S., 201 1, ApJ, 736, 2 
Ohsuga K., Mori M., Nakamoto T., Mineshige S., 2005, ApJ, 628, 

368 

Okuda T., Teresi V., Toscano E., Molteni D., 2005, MNRAS, 357, 
295 

Onken C. A., Ferrarese L., Merritt D., Peterson B. M., Pogge 
R. W., Vestergaard M., Wandel A., 2004, ApJ, 615, 645 

Peng C. Y., Impey C. D., Rix H.-W., Kochanek C. S., Keeton 
C. R., Falco E. E., Lehar J., McLeod B. A., 2006, ApJ, 649, 616 

Poindexter S., Morgan N., Kochanek C. S., 2008, ApJ, 673, 34 

Pooley D., Blackburne J. A., Rappaport S., Schechter P. L., 2007, 
ApJ, 661, 19 

Poutanen J., Lipunova G., Fabrika S., Butkevich A. G., Abol- 

masov P., 2007, MNRAS, 377, 1 187 
Proga D„ Kallman T. R., 2004, ApJ, 616, 688 
Proga D., Stone J. M., Kallman T. R., 2000, ApJ, 543, 686 
Revnivtsev M., Gilfanov M., Churazov E., Sunyaev R., 2002, 

A&A, 391, 1013 
Shakura N. I., 1972, AZh, 49, 921 
Shakura N. I., Sunyaev R. A., 1973, A&A, 24, 337 
Turnshek D. A., 1984, ApJ, 280, 51 
Vestergaard M., Peterson B. M., 2006, ApJ, 641, 689 
Vierdayanti K, Mineshige S., Ueda Y, 2010, PASJ, 62, 239 
Wambsganss J., 2006, in G. Meylan, P. Jetzer, P. North, P. Schnei- 
der, C. S. Kochanek, & J. Wambsganss ed., Saas-Fee Advanced 
Course 33: Gravitational Lensing: Strong, Weak and Micro Part 
4: Gravitational microlensing. pp 453-540 



© 2009 RAS, MNRAS OOO.mnTI 



